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Abstract 



We perform a high-statistics simulation of the three-dimensional randomly 
dilute Ising model on cubic lattices with L < 256. We choose a particular 
value of the density, x = 0.8, for which the leading scaling corrections are 
suppressed. We determine the critical exponents, obtaining u = 0.683(3), 
r] = 0.035(2), P = 0.3535(17), and a = —0.049(9), in agreement with previous 
numerical simulations. We also estimate numerically the fixed-point values 
of the four-point zero-momentum couplings that are used in field-theoretical 
fixed-dimension studies. Although these results somewhat differ from those 
obtained using perturbative field theory, the field-theoretical estimates of the 
critical exponents do not change significantly if the Monte Carlo result for 
the fixed point is used. Finally, we determine the six-point zero-momentum 
couplings, relevant for the small-magnetization expansion of the equation of 
state, and the invariant amplitude ratio R^ that expresses the universality of 
the free-energy density per correlation volume. We find R'^ = 0.2885(15). 

PACS Numbers: 75.10.Nr, 75.40.Cx, 75.40.Mg, 64.60.Fr 
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I. INTRODUCTION. 



During the last few decades many theoretical and experimental studies have investigated 
the critical properties of statistical systems in the presence of quenched disorder. Typi- 
cal examples are randomly dilute uniaxial antiferromagnets, for instance, Fe2:Zni_-EF2 and 
Mnj.Zni_j;F2, obtained by mixing a uniaxial antiferromagnct with short-range interactions 
with a nonmagnetic material. Experiments find that, for sufficiently low impurity concentra- 
tion 1-x, these systems undergo a second-order phase transition at Tc{x) < Tc{x = 1), with 
critical exponents independent of the impurity concentration. The experimental results have 
been summarized in Ref. [1], which reports a — —0.10(2), u — 0.69(1), and /3 = 0.350(9). 
These estimates are definitely different from the values of the critical exponents of the pure 
Ising universality class, where, e.g., a = 0.1096(5) (Ref. [2]), and thus indicate that the 
impurities change the nature of the transition that belongs to a new random universality 
class. In the presence of an external magnetic field, dilute uniaxial antiferromagnets show 
a different critical transition, belonging to the universality class of the random-field Ising 
model [3-6]. 

A simple model for dilute uniaxial systems is provided by the three-dimensional random 
Ising model (RIM) with Hamiltonian 



where the sum is extended over all nearest-neighbor sites, are Ising spin variables, and pi 
are uncorrelated quenched random variables, which are equal to one with probability x (the 
spin concentration) and zero with probability 1 — x (the impurity concentration) . For small 
1 — x, i.e. above the percolation threshold of the spins, this model shows a critical transition 
analogous to that observed in experiments and whose nature has been the object of many 
theoretical studies, see, e.g., Refs. [7-10]. 

Numerical Monte Carlo simulations [11-17] had long been inconclusive in setting the 
question of the critical behavior of the RIM. While the measured critical exponents were 
definitely different from the Ising ones, results apparently depended on the spin concentra- 
tion, in disagreement with renormalization-group (RG) theory. Only recently has the ques- 
tion been clarified. Ref. [13] showed the presence of very strong concentration-dependent 
scaling corrections with exponent uj — 0.37(6). Only if they are properly taken into account, 
the numerical estimates of the critical exponents become dilution independent as expected. 
Their final estimates are u = 0.6837(53) and r) = 0.0374(45), from which one also derives 
(3 = 0.3546(28) and a = —0.051(16) using scaling relations. These results arc in good agree- 
ment with the experimental ones reported above, although the numerical estimate of a is 
slightly different. 

Randomly dilute Ising systems can also be studied by using the field-theoretical (FT) 
approach [9,10]. The starting point is the cubic-symmetric Hamiltonian [18] 



where 0i is an iV-component field. By using the standard replica trick, it can be shown 
that in the formal limit N ^ such a model corresponds to a system with quenched 
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disorder effectively coupled to the energy density, as is the case of the RIM [18]. As is 
well known, the limit — ^ is a subtle one. In the standard pcrturbative approaches, 
the limit is taken naively — we simply set A^ = in the perturbative expansion — implicitly 
assuming that the replica symmetry is not broken. In recent years, however, this assumption 
has been questioned [19] on the ground that the RG approach may not take into account 
other local minimum configurations of the random Hamiltonian (1.2), which may cause the 
spontaneous breaking of the replica symmetry. However, a fixed-dimension perturbative two- 
loop calculation [20] in a perturbative approach proposed in Rcf. [19] finds that the standard 
replica-symmetric fixed point is stable under any replica-symmetry breaking perturbation, 
thereby supporting the standard approach. In this paper, we do not further consider this 
issue and in the following we always assume that the standard approach is correct. Note that 
the good agreement among numerical and field-theoretical results supports this assumption, 
although one cannot exclude that replica-symmetry breaking effects can only be seen very 
close to the critical point. 

In the FT approach one looks for stable fixed points in the region uq < (or, equivalently 
u < 0). If the pure fixed point at -u = is stable, disorder is irrelevant, while the presence 
of a new stable fixed point with u < 0, indicates that disorder is relevant and that dilute 
systems belong to a new universality class. Numerical and experimental results indicate 
that in dilute Ising systems the correct scenario is the second one and thus a new random 
fixed point should be present with u* < 0. 

The most precise FT estimates of critical quantities are presently obtained by using 
perturbative methods. However, in the case of random systems the perturbative approach 
faces new difficulties: the perturbative series are not only divergent, but are also non- 
Borel summable [21,22]. This means that even the knowledge of the complete perturbative 
series does not allow the exact computation of the critical quantities. These difficulties are 
clearly present in the ye-expansion and in the related minimal-subtraction scheme without e- 
expansion [23-25]. The expansion in ^ is not well behaved and does not allow quantitative 
determinations of the critical exponents, while in the minimal-subtraction scheme results 
are very sensitive to the resummation method. If the Chisholm-Borel method is used [24], 
no random fixed point is found with the longest available series (five loops). Apparently, 
four-loop series provide the most accurate results and increasing the length of the expansion 
does not help improving the precision of the results. On the other hand, if a double Pade- 
Borel resummation is used as proposed in Ref. [22], a random fixed point is found also at 
five loops [25]. The estimates of the critical exponents are in any case not very precise, 
and moreover, at variance with the fixed-dimension approach described below, the stability- 
matrix eigenvalues turn out to be complex. 

The fixed- dimension perturbative expansion in powers of two independent zero- 
momentum quartic couplings u and v (directly related to Uq and Vq defined in Eq. (1.2)) is 
apparently better behaved. Up to six loops, a random fixed point u*, v* is always found, 
although the estimates of u* and v* vary significantly with the order and the resummation 
method. In spite of that, the estimates of the critical exponents are quite precise, due to a 
relatively large insensitivity of the results to the position of the fixed point. The analysis of 
the six-loop series gives [26] u = 0.678(10), r] = 0.030(3), /S = 0.349(5), and a = -0.034(30). 
The agreement with the experimental and numerical results is again quite satisfactory; only 
the estimate of a seems to be slightly larger than the experimental result. 
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In this paper we present a new numerical study of the RIM. The purpose is to extend and 
possibly improve the numerical results of Ref. [13]. We estimate the critical exponents and, 
in particular, a in order to verify if the apparent small discrepancy between experiments 
and numerical results is really there. Moreover, we determine the four-point and the six- 
point zero- momentum couplings, and the universal ratio R^. As a byproduct we are able 
to check the accuracy of the FT approach by comparing Monte Carlo and FT estimates of 
the fixed-point values u* and v* . 

We have performed a high-precision Monte Carlo simulation of the model with Hamil- 
tonian (1.1) at J = 1 and density x — 0.8. Such a value has been chosen on the basis of 
the results of Ref. [13], where it was shown that scaling corrections are particularly small 
for such a value of x. This is fully confirmed by our analysis: We do not observe scal- 
ing corrections with exponent uj = 0.37(6), the correction-to-scaling exponent observed in 
Ref. [13] for generic values of the density. [27] Note that the absence of corrections with 
exponent uj also implies the absence of corrections with exponents 2a;, Sou, ... Therefore, 
we expect corrections to scaling with next-to-leading exponent 002 (<^2 = 0.8(2) according 
to field theory [28]). Unexpectedly, also these corrections are small. The RIM at x = 0.8 
is therefore an "improved" model [29-32,2], i.e. a model in which the leading correction to 
scaling is (approximately) absent in the expansion of any observable near the critical point. 

First of all, we determine the critical exponents by using two different methods. A 
first estimate is obtained by employing the extrapolation method of Refs. [33-35] (similar 
methods have been discussed in Refs. [36,37]). It allows us to determine the critical exponents 
from the high-temperature behavior of the susceptibility and of the correlation length. We 
also use direct finite-size scaling (FSS) methods, obtaining consistent estimates. Our final 
results are 

z/ = 0.683(3), (1.3) 
7] = 0.035(2), (1.4) 

from which, using scaling and hyperscaling relations, we obtain 

7 = i/(2- 77) = 1.342(6), (1.5) 

/3 = ^(1 + 77) = 0.3535(17), (1.6) 

5=^ = 4.80(11), (1.7) 

a = 2 - 3i/ = -0.049(9). (1.8) 

Our results are in good agreement with those of Ref. [13] and, in particular, confirm the 
discrepancy between the experimental and theoretical estimates of a. 

We also carefully check the validity of the hyperscaling relation 2 — a = 3z/. We analyze 
the specific heat at the critical point obtaining a/u — —0.115(28). Using the estimate (1.3) 
for Xjv we obtain 

2 n 

^ ~ ^ = 3.04(3), (1.9) 

which is fully consistent with 3. We also perform another check of hyperscaling, analyzing 
the specific heat and the energy at the critical point. We obtain 
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2 ni 

^-^ = 2.93(6), (1.10) 

again consistent with 3. 

Beside the critical exponents we also measure the four-point couplings G\ and G22 defined 
in Eqs. (A12) and (A14), which can be directly related to the fixed-point values u* and v*: 
Gl = V* and = u* We obtain: 

G\ = 43.3(2), 

G;, = -6.2(1) . (1.11) 

These estimates difi^er significantly from those reported in Ref. [26], which were obtained 
from the analysis of perturbative six-loop series: G^ = 38.0(1.5) and G22 = —4.5(6). Clearly, 
the non-Borel summability of the perturbative expansions gives rise to a large systematic 
error. It is also possible that the nonanalyticity of the RG functions [38-41] near the random 
fixed point plays an important role. 

These discrepancies on the estimates of u* and v* call for a reanalysis of the perturbative 
expansions of the critical exponents. By using the Monte Carlo estimate of u* and v* we find 
1/ — 0.686(4), Tj = 0.026(3), and 7 = 1.355(8). These estimates do not differ significantly 
from those obtained in Ref. [26] and are also in satisfactory agreement with the Monte 
Carlo results. Clearly exponents are quite insensitive to the exact location of the fixed 
point. We also try a different method for estimating critical quantities. It is based on an 
expansion around the Ising fixed point. Results are similar: u = 0.690(8), rj = 0.0345(20), 
and 7 = 1.355(10). Note that the estimate of f] is now in perfect agreement with the Monte 
Carlo result. 

In this paper we also determine some other universal amplitude ratios that involve high- 
temperature quantities. First, wc determine the six-point universal ratios rg, CI2, and C222, 
defined in Eq. (A13). The cocflicicnt rg is particularly important since it parametrizes the 
small-magnetization expansion of the critical equation of state in the high-temperature phase 
We find 



r* = 0.90(15). (1.12) 
Finally, we compute the universal ratio defined by 

i?+ ^ (aA+)i/V+, (1.13) 

where and /+ are defined in terms of the singular behavior of the specific heat C and of 
the correlation length { , Cgmg ~ ^ pa f^t'" for t = /?c — /9 — >■ 0+. We obtain 

= 0.2885(15), (1.14) 

in good agreement with other theoretical results [42]: R'^ = 0.290(10), obtained from the 
analysis of the corresponding six-loop perturbative series, and R^ = 0.282(3), derived from 
a quite precise approximation of the equation of state. 

The paper is organized as follows. In Sec. II we present the Monte Carlo results. In 
Sec. II B we determine the critical temperature by performing a careful analysis of the finite- 
size behavior of some RG invariant ratios near the critical point. In Sec. II C we determine 
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the four-point and six-point couplings by using the extrapolation method of Refs. [33-35]. In 
Sec. 11 D and 11 E wc determine v and ri by using again the extrapolation method and by also 
performing a more direct FSS analysis. Then, in Sec. 11 F we study the finite-size behavior 
of the energy and of the specific heat near the critical point. We obtain an independent 
estimate of a, that allows us to check the validity of the hyperscaling relation 2 — a — ?)V. 
Finally, in Sec. 11 G we compute the universal ratio K^. For this purpose, we generalize the 
extrapolation method of Refs. [33-35] to the energy. In spite of the necessary subtractions, 
the method works quite well, providing a rather precise estimate. In Sec. Ill we reanalyze 
the six-loop perturbative series of Ref. [26], using the new Monte Carlo estimate of the fixed 
point. We employ the different resummation methods discussed in Ref. [26] and also a new 
method based on an expansion around the Ising fixed point. Finally, in the appendix we 
report the definitions of the quantities that are used throughout the paper. 



II. NUMERICAL RESULTS 

A. The Monte Carlo simulation 

We have performed a high-precision Monte Carlo simulation of the model with Hamilto- 
nian (1.1) with J = 1 at density x = 0.8. Such a value has been chosen on the basis of the 
results of Ref. [13], who showed that for such a value of x scaling corrections arc particularly 
small. In the simulations we have considered cubic lattices of size L^, L = 16, 32, 64, 128, 
and 256, with periodic boundary conditions. Simulations have been performed for several 
values of j3 between 0.275 and 0.28578. Two thirds of the simulations refer to the interval 
0.275 < (3 < 0.2856 (we will call the corresponding data the high-temperature results), while 
one third of the CPU time was used in simulations in a narrow interval around the critical 
point, 0.28572 < j3 < 0.28578. The average number of samples for each (3 and L has been 
approximately 7 • 10^ (L = 16), 36 • 10=^ (L = 32), 27 • 10=^ (L = 64), 12 • 10=^ (L = 128), 
and 3 • 10^ (L = 256). The runs were performed on a cluster with Dual Athlon MP 1.2MHz 
processors. The total CPU time is approximately 17.4 CPU- years of a single processor. 
As random number generator we have used a combination of the Parisi-Rapuano generator 
[43] and of a congruential generator [44]. Results for each sample have been obtained as 
follows. Starting from a random spin configuration, we perform 2000 iterations, each of 
them consisting alternatively of a Metropolis sweep and of a full Swendsen-Wang update. 
Wc use both a local and a nonlocal dynamics to guarantee equilibration of short-distance 
and long-distance modes. Then, we perform 2000 full Swendsen-Wang updates, measur- 
ing all quantities (see appendix for definitions) every 10 iterations. To estimate correlation 
functions we use improved estimators that significantly reduce the statistical errors. Note 
that we have been much more conservative than Ref. [13]: There, only 200 Swendsen-Wang 
iterations were performed for equilibration. For quantities that involve the noise average 
of products of sample averages, there is a bias due to the finite length of the run for each 
sample. In order to take this bias into account we have performed a bias correction following 
Ref. [45]. 
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TABLE I. Results of the fit R{(5, L) 
freedom of the fit. 



R* + a{(5 - (5c)L^/''. DOF is the number of de grees of 



L 



'min 


x7dof 














16 


115.2/25 


1.6337(4) 


0.2857520(7) 


32 


17.5/20 


1.6385(6) 


0.2857477(6) 


64 


7.4/12 


1.6407(11) 


0.2857462(9) 






u& 




16 


134.4/25 


3.2348(22) 


0.2857530(8) 


32 


18.8/20 


3.2597(28) 


0.2857480(6) 


64 


7.4/12 


3.2710(54) 


0.2857463(9) 






U22 




16 


35.0/25 


0.1500(8) 


0.2857266(141) 


32 


25.0/20 


0.1484(6) 


0.2857425(68) 


64 


21.2/12 


0.1480(10) 


0.2857454(83) 










16 


22.6/25 


0.5921(2) 


0.2857414(6) 


32 


17.4/20 


0.5926(3) 


0.2857423(6) 


64 


11.8/12 


0.5934(6) 


0.2857432(9) 



0.804(44) 
0.727(43) 
0.726(60) 

0.817(46) 
0.727(43) 
0.726(60) 

1.23(51) 
1.06(54) 
0.93(66) 

0.733(34) 
0.708(37) 
0.722(54) 



B. Determination of the critical temperature 

As a first step in our analysis we have determined the critical temperature f3c- For this 
purpose we consider the results of the simulations for 0.28572 < {3 < 0.28578, which is 
a small interval around Pc- We consider four invariant ratios, t/4, U^, U22, and ^/L, see 
appendix for definitions. Standard FSS predicts that, in the FSS limit /3 — > /3c, — > 00 at 
{P — /3c)L^^'' fixed, each quantity i?(/3, L) behaves as 

R{(3,L) = R[{(3-(3,)L'/% (2.1) 

where R{z) is a universal function. Since (3 — (3c is particularly small for the data, we can 
expand R{z) in powers of z, keeping only the first term (we checked that the addition of the 
term of order z"^ does not change the results). Thus, we fit each quantity R{P, L) by using 

RiP,L)^R* + aiP-f3c)L'/'', (2.2) 

with R* , Pc, and u free parameters. In each fit we include all data with L > Lmin and, in 
order to detect corrections to scaling, we use L^i^ = 16, 32, and 64. The results are reported 
in Table I. 

There is a slight evidence of corrections to scaling, but it is interesting to note that they 
have opposite sign in C/4, Uq, and ^/L. Conservatively, we would obtain 

(3c = 0.2857447(24), (2.3) 



7 



where the central value is the average of the estimates obtained for U4 and ^/L (Lmm = 64), 
and the error is such to include one error bar for both U4 and C,/L. 

These fits are not particularly sensitive to the exponent z/, which is quite poorly deter- 
mined. We obtain 1/ — 0.72(6). One could imagine of improving the results by fixing i/. 
However, the dependence of (3c on 1/ is very small and no significant change is observed. 

In order to include scaling corrections we also perform fits of the form 

R{f3, L) = R* + a(/3 - (3^)1^^'' + bL'^, (2.4) 

where we include the leading scaling correction. These fits are not sensitive to the value of 
u and thus we fix it, taking u = 0.69. We keep a; as a free parameter, since we do not know 
which is the most important correction to scaling for our data. Indeed, the leading correction 
has exponent u; — 0.37(6), but there is evidence that for p — 0.8 leading corrections have a 
very small amplitude [13]. In order to be able to keep c<j as a free parameter, we analyzed at 
the same time two different observables. We restrict our attention to f/4, Uq, and ^/L, since 
U22 is too noisy. Using all data with L > 16 (if the data with L = 16 are discarded the fit is 
unstable) and taking properly into account the covariance between the two observables, we 
obtain: 

(a) analysis of ^/L and C/4: u = 0.70(11), (3c = 0.2857435(8), (^/L)* = 0.5943(8), = 
1.6502(24), b{^/L) = -0.017(7), b{U^) = -0.13(3); = 47.5, DOF= 50. 

(b) analysis of ^/L and Uq: uj = 0.71(11), (3c = 0.2857435(8), {^/L)* = 0.5942(8), = 
3.318(11), b{^/L) = -0.017(7), b{U4) = -0.67(16); = 46.7, DOF= 50. 

Here DOF is the number of degrees of freedom. The results are quite stable, indicating the 
presence of corrections with exponent u = 0.7(1), in agreement with the idea that scaling 
corrections with exponent lo ^ 0.4 are very small. The effective exponent to — 0.7(1) is 
quite close to the next-to-leading exponent predicted by perturbative field theory [28], i.e. 
UJ2 = 0.8(2). Thus, we mainly observe next-to-leading corrections, which in any case are 
quite small. In particular, they are of little relevance for C,/L. The coefficient b, cf. Eq. (2.4), 
is very small and the estimates of (3c obtained from the combined fits are fully compatible 
with those obtained for ^/L without scaling corrections. 

These analyses that keep into account scaling corrections hint at values of (3c lower than 
the estimate (2.3). We are thus led to consider 

(3c = 0.285744(2) (2.5) 

as our final estimate. 

From the above-reported analyses, we also obtain estimates of the invariant ratios R* at 
the critical point. We obtain 

0.5943(9), 

1.650(9), 
3.32(5), 

0.1480(10). (2.6) 
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We quote the results obtained in the fits (2.4), while the error is such to include also the 
result of the fit without scaling corrections and I/min = 64. For U22 we only consider the fits 
without scaling corrections. Note that U22 7^ 0, indicating the absence of self-averaging at 
the critical point, in agreement with the theoretical arguments of Ref. [46]. 

We can compare our results with those of Ref. [13]. They found Pc — 0.2857421(52), 
= 1.653(20), U^2 = 0.145(7), (^/L)* = 0.598(7), which are in full agreement with our 
final results. 



C. Determination of the four-point and six-point couplings 

In this Section we wish to determine the four-point couplings Gl and G22, and the 
six-point universal ratios rg, C42, and CI22, see appendix for definitions. Note that these 
quantities are defined in the high-temperature phase and one should take first the infinite- 
volume limit and then the limit j3 j3c, cf. Eq. (A14). In order to perform this task we 
have applied the extrapolation method of Rcfs. [33-35] to our high-temperature data, i.e. 
to the results with /3 < 0.2856 (corresponding to ioo{P) ^ 89). This method is extremely 
powerful in order to compute the infinite- volume behavior of critical quantities and it has 
been applied to several models, including spin glasses [47]. 

The idea is the following. Given a long-distance quantity S{(3,L), in the FSS limit we 
can write 

= Fs W,L)/L] + 0{L-,rn, (2.7) 

where s is an arbitrary (rational) number (here we always consider s = 2). Here Fs{z) is 
a universal function defined for < ^ < ^* = (^/L)* = 0.5943(9) such that Fs{0) = 1 and 
Fs{z*) = s°", where a is the exponent characterizing the critical behavior of S{f3,L) at the 
critical point, i.e., S{(3c, L) L'^ . Eq. (2.7) is the basis of the extrapolation technique since, 
in the absence of scahng corrections, it allows us to compute S'(/3, sL) on a lattice of size sL 
in terms of quantities defined on a lattice of size L and of the function Fs{z). In practice, one 
works as follows. First, one performs several runs, determining S{3,sL), S{f3,L), C,{l3,sL), 
and ^{j3, L). By means of a suitable interpolation, this provides the function F{z) for S and 
^. Then, S^{13) and ^qo{(3) are obtained from S{P, L) and ^(/3, L) by iterating Eq. (2.7) and 
the corresponding equation for ^{[3,L). Of course, one must be very careful about scaling 
corrections, discarding systematically lattices with small values of L till the results become 
independent of L within error bars. 

Let us first discuss the four-point couplings for which we will obtain quite precise es- 
timates. In Figs. 1, 2, and 3 we report the data for S{l3,2L)/ S{(3, L) for G4, G22, and ^, 
together with a fit of the data with L > 64. As discussed in Refs. [33,34], we parametrize 
Fs{z) with a polynomial [48] in e"^^^ of order n, increasing n until the of the fit changes 
by less than one by going from n to n + 1. In these analyses we have taken n = 7. The 
parametrization of Fs{z) as a polynomial in e~^/^ is theoretically motivated: indeed, for zero- 
momentum quantities, Fs{z) approaches 1 with corrections of order e~"/^, a ~ 1, as 2; — > 0. 
This choice is not strictly correct for ^ since in this case [49] F^{z) = 1 + 0(2;^). However, 
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0.1 0.2 0.3 0.4 0.5 



^(|3,L)/L 

FIG. 1. Ratios G4(/3, 2L)/G4(/?, L) vs ^{l3,L)/L. The solid curve is a fit using all data with 
L > 64. 



these power corrections are expected to be very small for our definition of finite- volume cor- 
relation length [50] , and therefore the systematic error due to our choice of parametrization 
should be small. 

Looking at the figures, it is quite difficult to distinguish any correction to scaling, i.e. 
systematic deviations from the fitted curve. However, at a closer look one may see that 
some points with L — 16 are out of the curve (in all cases by less than three error bars, so 
that these differences are barely significant). Conservatively, we have decided to discard all 
L — 16 data. In order to check further for corrections to scaling we have computed infinite- 
volume estimates S^{L) using only data with L > L^in, -^min = 32, 64. Additionally, we 
have also systematically discarded points that are far from the critical point by including 
only data with P > P^nin for several values of /?min- 

Using the extrapolation procedure we have outlined above, for each Lmin and /3min we 
obtain infinite- volume estimates ^ooiP), G4:,oo{P), and G22,oo(/3). To obtain the estimate at 
the critical point, the extrapolated values for the coupling constants have been fitted by 
using 

SM = S* + a{P,-P), (2.8) 

with l3c = 0.285744(2). The results are reported in Tables II and III for 0.275 < /^min < 0.284 
(corresponding to 4.45 ^ ^oo ^ 15.86). 

To check for corrections to scaling, we have also performed a different analysis. First we 
fit S{P, 2L)/S{P, L) taking into account a scaling correction with exponent cu, i.e. assuming 

= Fs m L)/L) + ^Gs m, L)/L) , (2.9) 

and use both functions Fs{z) and Gs{z) to perform the infinite- volume extrapolation. Then, 
the infinite-volume results for the coupling constants are fitted by using 
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TABLE II. Results for the renormalization constant G|. On the left we report the results for 
-^min = 32, on the right those for -Lmin = 64. We report two different x^- The first one (Xestr) 
refers to the fit that allows the determination of the curve Fc^iz), cf. Eq. (2.7), the second one 
(xlt) to the fit (2.8). Beside the we also report the number of degrees of freedom (DOF). The 
results have two errors: the first one is the statistical error, the second one gives the variation of 
the estimate as /?c is varied within one error bar, cf. Eq. (2.5). 



/5min 


xLtr/DOF 


Xit/DOF 


Gl 


Xitr/DOF 


Xit/DOF 


Gl 


0.2750 


14.9/22 


17.3/16 


43.65(6+0) 








0.2780 


14.9/22 


16.8/15 


43.63(6+0) 


9.0/11 


15.9/15 


43.56(9+2) 


0.2800 


14.9/21 


13.4/14 


43.59(7+0) 


9.0/11 


13.1/14 


43.51(10+0) 


0.2810 


14.9/20 


11.9/13 


43.60(8+1) 


9.0/11 


12.7/13 


43.52(10+1) 


0.2820 


14.8/19 


9.2/12 


43.55(8+1) 


9.0/11 


11.8/12 


43.47(11+0) 


0.2830 


12.1/17 


9.1/11 


43.41(11+1) 


7.4/10 


12.7/11 


43.45(12+0) 


0.2835 


12.0/15 


7.6/10 


43.34(12+1) 


7.4/9 


10.8/10 


43.40(15+1) 


0.2840 


9.1/13 


9.3/9 


43.34(17+2) 


7.1/8 


9.7/9 


43.32(18+1) 



TABLE III. Results for the renormalization constant G%n. Definitions are as in Table II. 



/3min 


xLr/DOF 


xit/DOF 




xLtr/DOF 


xit/DOF 




0.2750 


19.5/22 


13.6/16 


-6.18(3+0) 








0.2780 


19.5/22 


14.2/15 


-6.18(3+0) 


12.0/11 


11.3/15 


-6.25(4+0) 


0.2800 


19.1/21 


13.1/14 


-6.18(3+0) 


12.0/11 


11.5/14 


-6.25(4+0) 


0.2810 


17.4/20 


10.9/13 


-6.22(3+0) 


12.0/11 


11.3/13 


-6.25(5+1) 


0.2820 


17.4/19 


6.2/12 


-6.19(3+0) 


12.0/11 


8.9/12 


-6.23(5+0) 


0.2830 


16.9/17 


5.7/11 


-6.18(4+0) 


11.7/10 


8.4/11 


-6.23(6+0) 


0.2835 


16.8/15 


4.9/10 


-6.19(5+1) 


11.3/9 


8.2/10 


-6.25(6+1) 


0.2840 


12.1/13 


2.5/9 


-6.12(7+0) 


10.4/8 


1.5/9 


-6.11(8+1) 
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SM-S* + a{(5,-(3Y\ (2.10) 

In order to perform the analysis we should fix the exponents uj and v. We use v = 0.69, and 
repeat the analysis using cu — 0.8 (the next-to-leading exponent predicted by field theory 
[28]) and u — 0.4, which is the leading exponent determined in Ref. [13]. The results 
corresponding to c<j = 0.8 are reported in Tabic IV. It is essential to include the results 
with L = 16 in the analysis, otherwise the data do not show FSS corrections and the fit is 
unstable. Therefore, we cannot check the goodness of the Ansatz by discarding data with 
small L, i.e. present results for different values of Lmin as done before. 

Let us first discuss the results for G^. The fits without corrections to scaling show a 
significant (at the level of the reported errors) decrease as /^min is increased and also a slight 
dependence on Lmin- Corrections to scaling are positive and the estimate decreases with 
increasing Pram, so that one only obtains an upper bound Gl < 43.4. On the other hand, 
the fit with cu — 0.80 gives results independent of /^min within error bars; moreover the 
of the fit (2.10) is systematically lower than that of the fit (2.8). Clearly the data are very 
well fitted by assuming a correction-to-scaling exponent u = 0.80. For uj = 0.40 the results 
strongly depend on (3,^i^, varying from Gl = 42.49(14) for P^nm = 0.275 to Gl = 43.47(40) 
for /3inin = 0.2835. Also, the is larger than the obtained using uj = 0.8. There is 
therefore little evidence for such a small correction-to-scaling exponent, confirming again 
that for density x—0.8 the leading scaling corrections are very small. As final estimate we 
take 

= 43.3(2), (2.11) 

which is consistent with all results. 

Let us now consider the results for 6^22- The results of Table 111 show no dependence on 
Pmin and a tiny dependence on Lmin which could be of purely statistical origin. With the 
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0.1 0.2 0.3 0.4 0.5 

?(P,L)/L 

FIG. 3. Ratios 2L) L) vs ^(/?, L) /L. The solid curve is a fit using all data with L > 64. 

present error bars there is no evidence for nonanalytic scahng corrections and indeed the 
results obtained using u; — 0.80 (see Table IV) are perfectly consistent with those of the fits 
with purely analytic corrections. Our final estimate is 

= -6.2(1). (2.12) 

The error is rather conservative and is such to include all estimates. 

We should note that our results are compatible with the FT predictions of Ref. [28], where 
it is shown that in infinite volume G22,oo{P) and G4,oo(/9) have scaling corrections with next- 
to-leading exponent lo — 0.8(2) of similar relative size. In particular, the results of Ref. [28] 
give aG22 = ~0.23(10)aG4 where a is the coefficient defined in Eq. (2.10). From the fits we 
obtain instead a^^ = 20(5) and 0^22 = —1.5(1.0). The errors should be taken with caution: 
They simply give the variation of the parameter a with /9min for 0.275 < /^mm < 0.281 (for 
larger values the statistical error becomes larger than \a\) and do not include any possible 
systematic effect. Assuming these values with their errors, we estimate a^ja = — 0.08(5)aG4, 
which is in reasonable agreement with the FT result. 

As a check we have repeated the analysis for H4 = G4 + 3G22- Since the procedure is 
nonlinear, this represents an important consistency check. We obtain 

7/4* = 24.7(2), (2.13) 

which is in full agreement with the estimates reported above. 

Finally, we consider the six-point couplings Tq, C42, and C222- We apply again the 
extrapolation procedure we have used for G4 and ^22- However, in this case there are larger 
systematic errors. The extrapolation curve F{z), cf. Eq. (2.7), is poorly determined for 
z ^ 0.3, since the six-point couplings have large statistical errors when ^{f3,L)/L is small. 
A large error on the curve gives a large systematic error on the extrapolations and induces 
correlations among the results for different /3 (such correlations are instead small for G4 



13 



TABLE IV. Results for G| and G22 for fits with a correction-to-scaling exponent uj = 0.8. We 
report two different x^- The first one (Xestr) refers to the fit that allows the determination of the 
curve F{z), cf. Eq. (2.9), the second one (xit) to the fit (2.10). Beside the we also report the 
number of degrees of freedom (DOF). We only report the statistical error; the error due to Pc is 
negligible. 



/3min 


xLtr/DOF 


Xit/DOF 


^4 


xLtr/DOF 


Xit/DOF 


•-^22 


0.2750 


30.0/29 


12.6/16 


43.28(8) 


31.8/29 


12.1/16 


-6.15(3) 


0.2780 


29.9/28 


10.9/15 


43.31(9) 


31.6/28 


11.1/15 


-6.16(4) 


0.2800 


29.5/26 


10.8/14 


43.30(10) 


31.1/26 


9.9/14 


-6.16(4) 


0.2810 


25.4/24 


9.0/13 


43.43(11) 


22.1/24 


8.8/13 


-6.26(5) 


0.2820 


24.3/22 


9.0/12 


43.30(13) 


22.0/22 


5.8/12 


-6.21(6) 


0.2830 


21.3/19 


6.4/11 


43.37(17) 


21.1/19 


3.6/11 


-6.25(7) 


0.2835 


13.4/16 


7.1/10 


43.53(22) 


18.1/16 


3.2/10 


-6.33(9) 


0.2840 


10.6/13 


7.9/9 


43.26(29) 


12.3/13 


2.7/9 


-6.00(13) 



and ^22)- Because of them, it is difficult to set correct error bars and to determine the 
final results. For these reasons we have taken a conservative attitude. We have generated 

different sets of extrapolated data by changing the degree n of the interpolation polynomial 
and /9inin (in all cases we set Lmin = 16 in order to have a sufficiently large number of points 
in the small- 2; region). Then, we determine the range that includes most (approximately 
2/3) of the extrapolated data with their error bars. The central value of such an interval 
gives the final result, while its half width gives the error. We obtain 

= 0.90(15), 
C:2 = 0.12(5), 

C*22 = 0.45(15). (2.14) 

In Fig. 4 we report the final results together with the results for a single extrapolation, so 
that the reader can see the quality of our numerical results and how much conservative our 
final error bars are. Note that in this analysis we have not made any attempt to evaluate 
corrections to scaling. In any case, we expect them to be small compared with the large 
errors we quote. 

The numerical estimates can be compared with the FT results. We shall discuss the 
four-point couplings in Sec. III. For what concerns the six-point couplings, the analysis of 
the available four-loop series in the fixed-dimension scheme [51] is not very precise. We only 
mention the estimate — l.lloi reported in Ref. [42], which agrees with the more precise 
result (2.14). 

D. The exponent u 

In order to determine the exponent u we have followed two different strategies: (a) we 
have determined extrapolated values ^00 (Z^) using our extrapolation method and then we 
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FIG. 4. Infinite- volume results for the six-point ratios r^^P), C42(/3), and C222(/3) from the 
extrapolation with Lmin = 16, /3min = 0.275, and degree of the interpolation polynomial n = 7. On 
the horizontal axis we report ^oo{P)- We also report the final results Vq = 0.90(15), CI2 = 0.12(5), 
and CI22 = 0.45(15). 



have performed a fit ^oo{f^) ~ {Pc — P)^"', (b) we have performed a purely FSS analysis 
without extrapolations. The advantage of method (b) is that we can keep both /3c and a; as 
free parameters, and thus check the previously determined value for a-nd our claim that 
scaling corrections with small values of u are tiny. 

In the first case we proceed as before. If we neglect scaling corrections, we can generate 
infinite-volume estimates ^ooiP) by using the extrapolation precedure based on Eq. (2.7) and 
then we can determine z/ from fits of the form 

In Coo = -i^HPc -(3) + a + b{(3, - (5). (2.15) 

As a second possibility we can include corrections with exponent u in the extrapolation, see 
Eq. (2.9), and fit the extrapolated data with 

iniM = -^H(^c - /5) + a + h{(5, - pyr (2.16) 

We report in Fig. 5 the results for the analytic fit with Lmin = 32 and Lmin = 64 and for 
the nonanalytic fit with u — 0.80 and Lmin = 16 for several values of /5min- We observe that 
results with and without nonanalytic corrections differ significantly, much more than the 
statistical errors. However, corrections have opposite trends. Analytic fits give results that 
decrease as /3min is increased while nonanalytic fits have the opposite behavior. Compatible 
results are obtained for /3min > 0.283. For (3 = 0.283 we obtain v = 0.6847(15) from the 
first analysis with Lmin = 64 and p = 0.6818(24) from the second analysis. A reasonable 
estimate would be 0.683(3) which includes all results. It should be noted that the results 
depend strongly on (3^ that has been fixed to f3c = 0.285744(2). By varying (3^ within one 
error bar, u varies approximately by ±0.003. Thus, collecting everything together, this 
method gives the final result 
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FIG. 5. Results for the exponent v vs /3min- We report the results obtained by fitting the 
extrapolated high-temperature data without scaling corrections (they are labelled Lmin = 32 and 
Lmin = 64) and with scaling corrections with exponent u) = 0.8 (labelled oj = 0.80), and the results 
obtained by using all data and the parametrization (2.20) (labelled "FSS"). The horizontal lines 
correspond to the final result v = 0.683(3). 



1/ = 0.683(3 + 3). (2.17) 

The estimate (2.17) has been obtained by using only the high-temperature data, i.e. those 
with j3 < 0.2856. A more precise estimate can be obtained by performing a direct FSS 
analysis that allows us to use both the high-temperature and the critical-point results for 
the correlation length. In this way, we do not need to fix [3c and we can also keep a; as a 
free parameter. We start from the general FSS expression 

^^F[u^{t)L"\u2{t)L-% (2.18) 

where only one correction-to-scaling operator has been taken into account. Here t = {f3c—(3), 
and Ui{t) and U2{t) are nonlinear scaling fields satisfying Ui{0) = and 1*2(0) 7^ 0. Equation 
(2.18) can be rewritten in the scaling limit as 

I = MtL'/n + L-^'"" hitL^'n + L-^hitL^'"), (2.19) 

with /i(0) 7^ and /3(0) 7^ 0. The three scaling functions represent the three types of 
contributions we expect: fi{x) is the leading FSS curve, f2{x) corresponds to the analytic 
corrections, and /3(.t) is the nonanalytic FSS correction. The function f2{x) is related to 
fi{x) by f2{x) ~ f[{x)^ a relation that follows from the fact that this correction is due 
to the expansion of the scaling field Ui{t) = t + at^ + 0{t^). The existence of the infinite- 
volume limit fixes the behavior of these functions for a; — > 00: /i(a;) ~ x"^, f2{x) ~ x^~" 
and h{x) ~ x^'^-'^>. 
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TABLE V. Results of the FSS analysis using the parametrization (2.20) with n = 12 as a 
function of /3min- Here DOF is the number of degrees of freedom of the fit. 



xVdof 



UJ 



0.2750 
0.2780 
0.2800 
0.2810 
0.2820 
0.2830 
0.2835 



51/62 
50/60 
45/57 

43/54 
41/51 
39/47 
30/43 



0.6851(15) 
0.6842(22) 
0.6823(36) 
0.6803(52) 
0.6827(36) 
0.6836(28) 
0.6836(23) 



1.71(18) 
1.60(28) 
1.32(35) 
1.10(34) 
1.51(66) 
2.1(1.3) 
3.6(2.5) 



0.2857434(7) 
0.2857434(7) 
0.2857433(7) 

0.2857435(8) 
0.2857433(7) 
0.2857430(7) 
0.2857428(8) 



Now, we wish to determine f3c, i^, and uj by fitting the numerical data for ^/L with 
Eq. (2.19). For this purpose, we must parametrize the scaling functions reported above with 
expressions containing a finite number of parameters. Thus, we consider an even integer n 
and parametrize: 

I = [P„(tLV-) + L-(l +piLV-)--g„(iLV-)]"'/" , 

n n/2 

P„(x) = Yl Qn{x) = J2 hx^\ (2-20) 

i=0 i=Q 

where t = {j3c — (3). This parametrization has the correct behavior both for small and large 
X for any n. Note that our choice of writing Qn( polynomial in is due only to 

practical considerations: since corrections are small, Qn{x) cannot be parametrized with as 
many parameters as P„(x); on the other hand, to guarantee the correct large-x behavior 
Qn{x) and Pn{x) must have the same degree. We also made some analyses writing Qn{x) 
as a polynomial in — in this case n must be a multiple of 3 — without finding significant 
differences. The value of n has been chosen on the basis of the of the fit: n has been 
increased until the does not change significantly as n is increased by 2. The results of 
the analyses we present correspond to n = 12. 

In the analyses we have kept z/, cj, p, /9c, {oj}, and {hi\ as free parameters. The fit is 
stable only if we include the data with L = 16, otherwise there are no scaling corrections 
at the level of our error bars and therefore it is impossible to determine ui. To check for 
corrections, we have systematically discarded points far from P^, including in the analysis 
only data with (3 > Pram for several values of Pmin- In order to have a mostly linear fit, we 
have fitted (f )""/''. The results arc reported in Table V. 

The results for Pc are quite stable, in full agreement with the analyses reported in Sec. 
JIB for ^/L and with the final estimate (2.5). The estimates of cu vary significantly and 
have a quite large error. There is however very little evidence of scaling corrections with 
a; ^ 1, as already discussed in Sec. II B. Finally, let us consider u. As it can be seen from 
Fig. 5 the results are in rough agreement with those found before and apparently become 
independent of /Smm for P ^ 0.282. Our final estimate is 

u = 0.683(3), (2.21) 
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FIG. 6. Ratios ({p, 2L) /C{P, L) vs ^(/?, L) /L. The solid curve is a fit using all data with L > 64. 

which is compatible with the results of all analyses. The estimate (2.21) is in very good 
agreement with the result u — 0.6837(53) of Ref. [13]. 

E. The exponent rj 

In order to compute the exponent r] and correspondingly 7 = (2 — 77)2/, one can analyze 
the critical behavior of the susceptibility x- However, we have found much more convenient 
to analyze ( = xlC- Indeed, because of statistical correlations between x ^-nd ^^e relative 
error on C is significantly smaller than that on x- Moreover, the FT analysis of Ref. [28] 
indicates that C, has much smaller scaling corrections than x- 

In order to determine 77 we perform two different analyses. The first one uses the 
high-temperature data and follows closely the analysis of the correlation length presented 
in Sec. II D. Given CiP^L) we determine the curve F(^{z), cf. Eq. (2.7). The ratios 
C(/3, 2L)/C(/3, L) are reported in Fig. 6, together with a fit of the points with L > 64. 
There are clear scaling corrections, especially when ^(/3, L)/L ^ 0.3: the data with L — 16, 
and also some points with L = 32, are systematically above the curve. As before, we dis- 
card all points with L = 16 and perform the extrapolation using data with L > Lmiru with 
-^^min = 32, 64, in order to detect scaling corrections. The extrapolated values are fitted by 
using 

lnCoo(/9) = V^HPc - /?) + a + K/?c - /?) (2.22) 

and fixing Pc — 0.285744(2). The results are reported in Table VI. They show a systematic 
dependence on Lmin and /?min and seem to indicate rji/ ~ 0.0240, but it is difficult to set an 
error bar. Conservatively, we write 

771/ = 0.0240(15), (2.23) 
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TABLE VI. Results for the critical exponent riv. On the left we report the results for Lmin = 32, 
on the right those for Lniin = 64. We report two different x^- The first one (Xcstr) refers to the 
fit that allows the determination of the curve F^{z), cf. Eq. (2.7), the second one (xit) the fit 
(2.22). Beside the we also report the number of degrees of freedom (DOF). The results have 
two errors: the first one is the statistical error, the second one gives the variation of the estimate 
as Pc is varied within one error bar, cf. Eq. (2.5). 



/3min 


xLtr/DOF 






Xe'str/DOF 


Xit/DOF 


rju 


0.2750 


37.1/22 


16.2/15 


0.0249(1+1) 








0.2780 


37.1/22 


16.1/14 


0.0249(2+1) 


15.7/11 


12.4/14 


0.0253(2+1) 


0.2800 


36.7/21 


11.0/13 


0.0246(2+1) 


15.7/11 


11.6/13 


0.0252(2+1) 


0.2810 


36.7/20 


5.9/12 


0.0243(2+1) 


15.7/11 


10.6/12 


0.0251(3+1) 


0.2820 


36.3/19 


5.5/11 


0.0242(3+1) 


15.7/11 


11.0/11 


0.0250(3+1) 


0.2830 


35.1/17 


4.0/10 


0.0241(4+1) 


15.5/10 


9.9/10 


0.0248(4+2) 


0.2835 


33.2/15 


1.5/ 9 


0.0236(5+2) 


15.5/9 


4.9/9 


0.0243(6+2) 


0.2840 


22.6/13 


4.7/ 8 


0.0236(7+2) 


12.0/8 


2.7/8 


0.0245(8+1) 



that includes the estimates with /^min = 0.284 with their error bars. Using the estimate of v 
reported in Sec. II D, u — 0.683(3), we obtain 

ri = 0.035(2). (2.24) 

We have also repeated the analysis including scaling corrections with uj = 0.8. In the 
extrapolation procedure we use Eq. (2.9), and then we fit the extrapolated data by using 

lnCoo(/3) = V^HPc -f3) + a + 6(/3e - PT". (2.25) 

The results are reported in Table VII. They show a systematic upward trend, indicating 
that scaling corrections with u; — 0.8 do not correctly describe the data. This is also evident 
from the of the fit (2.25) which is quite large for small /3min- A lower value of a;, say 
ou — 0.4, gives even worse results, confirming that corrections with a; ^ 1 are very small, in 
agreement with the FT analysis of Ref. [28]. In any case the results of Table VII indicate 
that 7]iy > 0.0220, in agreement with Eq. (2.23). 

The second analysis we consider uses the data at the critical point. For /3 /3c we have 

C(/?, L) = L-^'MtL'/^) + L-^--h{tL''^), (2.26) 

where t = (3c — If we expand for i we can write 

In C(/9, L) = -r]\nL + a + htL^/'' + cL"'". (2.27) 

We first perform an analysis neglecting scahng corrections (c = 0) for several values of -Lmin- 
If we fix /5c = 0.285744(2) and v = 0.683(3) we obtain: 

Lrain = 16: 77 = 0.0331(2 + 6), xVDOF = 76/26; 
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TABLE VII. Results for the critical exponent ijiy using scaling corrections with exponent 
u; = 0.8. Wc report two different x^- The first one (Xcstr) I'cfcrs to the fit that allows the de- 
termination of the curve F(^{z), cf. Eq. (2.9), the second one (xit) ^^^^ (2.25). Beside the 
we also report the number of degrees of freedom (DOF). We only report two errors: the first one 
is the statistical error, the second one gives the variation of the estimate as Pc is varied within one 
error bar, cf. Eq. (2.5). 



/3min 


xLtr/DOF 


Xit/DOF 


rji' 


0.2750 


24.5/29 


174/15 


0.0179(2+1) 


0.2780 


24.5/28 


55.8/14 


0.0193(2+1) 


0.2800 


24.1/26 


29.0/13 


0.0203(3+1) 


0.2810 


22.5/24 


19.3/12 


0.0210(4+2) 


0.2820 


20.0/22 


11.0/11 


0.0217(5+2) 


0.2830 


19.4/19 


6.1/10 


0.0228(7+2) 


0.2835 


19.2/16 


4.8/9 


0.0226(8+2) 


0.2840 


14.2/13 


5.3/8 


0.0230(12+3) 



i^min = 32: = 0.0347(3 + 8), xVDOF = 27/21; 

^min = 64: = 0.0354(6 + 14), xVDOF = 22/13. 

The first error is statistical, while the second one is due to the error on /3c. The error due to 
the error on u can be neglected. There are apparently some scaling corrections, but it should 
be noted that the difference between the results with Lmin = 32 and 64 is not significant 
given the statistical error bars. We also perform an analysis with scahng corrections with 
exponent lo — 0.80. Using all data, i.e. taking Lmin = 16, we obtain 77 = 0.0382(8 + 16), 
c = -0.111(16 + 23), xVDOF = 29/25. This result is higher than the estimates obtained 
before, but still compatible with the quoted errors. 

The analysis at the critical point gives therefore results that are in full agreement with 
those obtained before, confirming the correctness of the estimate (2.24), and with the result 
of Ref. [13], 77 = 0.0374(45). 

F. The exponent a and hyperscaling 

In this Section we wish to determine the exponent a and check the hyperscaling relation 
2 — a = du. Unfortunately, we have measured the specific heat C(/3, L) only near the critical 
point and thus we can determine a only from the behavior for jS ^ [3c. 

RG predicts that, in the FSS hmit, the energy scales as 

E{(5, L) ^ Sbuik(/9) + L("-i)/-^i(tLV-) + L^^-^y--^g^{tL^'^), (2.28) 

where t = Pc — P and we have included one scahng correction. Note that £'buik(/3) is 
expected to have an exponentially small dependence on L which can be neglected for all 
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practical purposes. Near the critical point we can expand this expression in powers of tL^/'^ 
obtaining for the energy and the specific heat the expressions 

E{P, L)^aE + bEt + ceL^''-^^/'' + dEtL'^''' + egL^^-^)/'^-'^ , (2.29) 
C(/3, L)^ac + bet + ccV^''' + dctL^'^^'-^l^ + ecL"/""'", (2.30) 

with ac — bE and cc = dE- 

We wish now to determine ajv. We consider the specific heat, since in this case the 
singular behavior is stronger. In Eq. (2.30) we have not used hyperscaling, so that the 
expression depends on two independent exponents a and z/, or more precisely ajv and \ jv. 
In order to simplify the analysis, we now analyze the specific heat using hyperscaling to 
rewrite the correction term. Thus, Eq. (2.30) becomes 

C(/3, L) = ac + hct + ccL"/'^ + dctL'/''^''/''^'^ + e^L"/'^-'^, (2.31) 

that depends only on ajv. We begin by neglecting the correction-to-scaling term, i.e. we 
set 6(7 — 0. Fixing (3c = 0.285744(2), a five-parameter fit of the data with L > Lmin gives: 

L^in = 16, a/i/ = -0.119(7 + 6), xVDOF = 12.0/24; 
L^in = 32, a/u = -0.115(17 + 11), xVDOF = 11.6/19. 

The first error is the statistical one, while the second one gives the variation of the estimate as 
(3c is varied by one error bar. Two things should be noticed: the analytic term proportional 
to t is very small compared with the statistical errors. Indeed, be = 23(83 -|- 8) and be = 
—30(167 -|- 24) in the two fits. Moreover, corrections to scaling are apparently small, since 
the results do not depend on Lmin. As a check, we have also performed an analysis with a 
correction term (ec 7^ 0). We find a/iy = -0.109(50 + 6), ec = -0.3(1.4), and xVDOF 
= 11.9/23, for u — 0.8 and Lmm = 16. The coefficient ec is compatible with zero, while 
the of the fit changes only by 0.1 in spite of the fact that we have added an additional 
parameter. There is no evidence of scaling corrections. As final estimate oi a/u we take the 
result with L^m = 32 and ec = 0, 

^ = -0.115(28). (2.32) 

If we use hyperscaling and the estimate of f of Sec. IID, 1/ — 0.683(3), we obtain a/u — 
—0.072(13). This estimate is in reasonable agreement with (2.32) confirming the validity 
of hyperscaling. Using Eq. (2.32) and u = 0.683(3) we obtain a — —0.079(19), which is in 
reasonable agreement with the estimate obtained using hyperscaling a = —0.049(9). 

We can also check hyperscaling directly by comparing the results for the energy and the 
specific heat. If we define 9i = a/v, 62 = {2 — a)/u, we can rewrite Eqs. (2.29) and (2.30) 
as 

E{(3, L)=aE + bEt + CijL^^i-^^)/' + dEtL'\ 

C{(3, L)=bE + dEL'^' + dctL^^'^'+^'^Z^, (2.33) 

where we have neglected scaling corrections and we have set 6^ = in agreement with the 
previous analysis. If hyperscaling is satisfied we should find 6*2 = 3. A combined analysis of 
E{(3, L) and C{(3, L) fixing (3^ gives 
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FIG. 7. Ratios R{P,2L)/R{P,L) vs ({P,L)/L. The solid curve is a fit using all data with 
L > 64. 

Linin = 16: a/i/ = -0.119(6 + 6), ^2 = 2.952(17 + 8), xVDOF = 33.4/51; 

Lmin = 32: a/u = -0.112(16 + 9), 02 = 2.930(43 + 12), xVDOF = 29.3/41. 

Here, as before, the first error is the statistical one, while the second gives the error due to 
(3c- The estimates of a/u are compatible with Eq. (2.32), while 62 is in reasonable agreement 
with 3, confirming again the validity of hyperscaling. 



G. The universal ratio 

We now wish to compute the universal ratio that is related to the universality of the 
singular part of the free energy per correlation volume. More precisely, let us consider the 
infinite- volume free energy density J-'oo{P) 

^00 = -^EZ, (2-34) 

and the infinite- volume specific heat Coo(/3) = 9^^oo(/3)/9/3^. In the critical limit t = 
/3c — /? ^ we can write 

^oo(/3) = J^hMP) + F+f-"" + corrections, 

C^{(3) = Cbuik(/?) + A+r'^ + corrections, (2.35) 

where = (2 - a)(l - a)F+. If ^ooiP) ~ f^t'" in the same limit t 0+, RG predicts the 
universality of F+(/+)'\ or equivalently of 

R^=(aA^y^'f+. (2.36) 
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TABLE VIII. Results for the universal ratio R'^. Definitions are as in Table II. The error due 
to Pc is negligible. 



/3min 




Xlt/DOF 




xLr/DOF 


Xit/DOF 




f) 2750 


14 7/1Q 


40 2/1fi 


28674('1 S^l 








0.2780 


14.7/19 


31.3/15 


0.28688(14) 


10.4/8 


20.1/15 


0.28677(16) 


0.2800 


14.6/18 


23.9/14 


0.28703(15) 


10.4/8 


16.5/14 


0.28689(17) 


0.2810 


14.6/17 


20.9/13 


0.28712(17) 


10.4/8 


15.8/13 


0.28698(19) 


0.2820 


13.8/16 


20.4/12 


0.28718(18) 


10.4/8 


13.2/12 


0.28709(20) 


0.2830 


13.6/14 


15.4/11 


0.28744(22) 


10.3/7 


10.8/11 


0.28726(23) 


0.2835 


12.9/12 


12.8/10 


0.28749(26) 


10.3/6 


9.45/10 


0.28736(27) 


0.2840 


12.9/10 


11.5/9 


0.28737(32) 


10.2/5 


8.63/9 


0.28730(33) 



We now compute using our high-temperature results for the energy E{(3, L) and for the 
correlation length ^{(3, L). We define a quantity R{I3, L), 



a 



■{E{p,L)-aE-bEt)t 



1/3 



mL), (2.37) 



1 — a 

where qe and Be are defined in terms of the expansion of ^buik(/5) for i — > 0: 



d(5 

It is easy to check that 



= £^buik(/9) = + hEt + 0{t'). (2.38) 



lim lim i?(/3, L) = Rj . (2.39) 



In order to compute R{P,L), we must specify the values of the two constants ge and 6e 
in Eq. (2.37). For this purpose we exploit the fact that -Ebuik(/3) is the same function 
in Eq. (2.28) and in Eq. (2.38), so that and Be coincide with the constants defined in 
Eq. (2.29). Thus, aE and Be can be determined independently by using the critical-point data 
for the energy and the specific heat. We thus proceed as follows. We consider Eq. (2.33), 
fix /3c = 0.285744(2), 02 = 3, a/iy = 2/iy - 3, ly = 0.683(3), and compute and Be by 
analyzing E{(3,L) and C{P,L) near the critical point. Then, we determine R{P,L). The 
error on R{I3,L) takes into account the error on E{(3,L), ^(/9,L), aE, and Be, and also the 
variation of the estimates as u and j3c vary within one error bar. In order to be conservative, 
we use a worst-error estimate summing all errors together. Once R{(3, L) has been computed 
we use the extrapolation method presented in Sec. II C. 

In Fig. 7 we report the ratios R{P,2L)/ R{/3, L) together with a fit of the data with 
L > 64 (a good fit is obtained by using a polynomial with n — 10). Apparently, there are 
no scaling corrections, but at a closer look one finds systematic deviations for L — 16. As 
before, these points will be discarded in the analytic fits. 

The results of the fits with analytic corrections are reported in Table VIII for Lmin = 32 
and Lmin = 64. For small /^min they show an upward trend and then apparently stabilize 
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TABLE IX. Results for the universal ratio R'^ using scaling corrections with exponent lo = 0.8. 
Definitions are as in Table IV. The error related to Pc is negligible. 



/3min 


Xestr/DOF 


Xfit/DOF 


Rt 


2750 


"io /25 


10 5/1fi 


28908 flQ") 


0.2780 


30.8/24 


9.5/15 


0.28897(21) 


0.2800 


30.5/22 


8.6/14 


0.28895(24) 


0.2810 


29.7/20 


7.6/13 


0.28897(26) 


0.2820 


28.5/18 


8.5/12 


0.28898(30) 


0.2830 


27.9/15 


8.5/11 


0.28905(36) 


0.2835 


17.8/12 


5.2/10 


0.28951(45) 


0.2840 


15.7/9 


6.4/9 


0.28967(56) 



around 0.2874(3). In order to check the role of the corrections to scahng, we have repeated 
the analysis by adding scaling corrections with exponent uj = 0.80. The results are presented 
in Table IX. They arc now substantially independent of (3^i^ confirming that the data are 
very well fitted by assuming such an exponent. The final estimate is somewhat higher than 
that obtained in the analytic fits, indicating that in this case nonanalytic scaling correction 
may play an important role. We do not know which of the two fits is the most reliable one 
and thus we have taken as final estimate 

i?J = 0.2885(15), (2.40) 

which is compatible with the results of both analyses. 

The estimate (2.40) is in good agreement with the results of other methods. A six- 
loop computation in the fixed-dimension FT approach gives [42] — 0.290(10), while, 
by using approximate parametric representations of the equation of state, one obtains [42] 
R+ = 0.282(3). 

III. COMPARISON WITH FIELD-THEORY RESULTS 

The critical behavior of the RIM has been extensively studied using the FT approach. 
Quantitative predictions can be obtained by using different techniques: perturbative meth- 
ods in the four-point renormalized couplings in fixed dimension d = 3 or in ^/e, e = 4 — d, or 
nonperturbative methods based on approximate RG equations, see Rcfs. [23-26,28,51-56] 
The most accurate results have been obtained in the first approach: six-loop expansions for 
the /3-functions and the critical exponents have been derived and analyzed in Refs. [57,26]. 
The corresponding estimates of the critical exponents, e.g., u — 0.678(10) and r) — 0.030(3), 
are in satisfactory agreement with the Monte Carlo results presented before. 

The main problem of the perturbative approach is the non-Borel summability of the series 
[21,22]. This fact makes the analysis more subtle and less precise than in the case of the pure 
Ising model. The difficulties of the perturbative approach appear in the determinations of 
the fixed-point values u* and v* of the renormalized couplings (they are normalized so that at 
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tree level u = Uo/m, v — Vo/m, m being the renormalized mass), which are directly related to 
the quantities we have measured in the Monte Carlo simulation (see Ref. [42] for a derivation 
of these relations): Gl = v*, G22 = '"*/3, and = Gl + 3G22 = u* + v*. The analysis 
of the six-loop series provided results somewhat dependent on the resummation method 
[58]. Indeed, we found [26]: Gl = 37.7(2), G*^^ = -4.3(6), = 24.8(1.8) (double Pade- 
Borel method); Gl = 36.8(3), = -4-0(1), = 24.8(6) (conformal Pade-Borcl method); 
Gl = 38.6(7), G*22 = -4.8(2), = 24.2(9) (direct conformal method). The estimates of 
are in good agreement with the Monte Carlo result (2.13), = 24.7(2). On the other hand, 
the estimates of Gl and G22 — combining the results we would have guessed Gl = 38.0(1.5) 
and G22 — —4.5(6) with errors that are, at first sight, quite conservative — differ significantly 
from the Monte Carlo estimates (2.11) and (2.12). 

These discrepancies call for a reanalysis of the perturbative series of the exponents, 
verifying if the use of the Monte Carlo results for Gl and G22 leads to significantly different 
estimates. We have thus repeated the analysis, using the different resummation methods 
outlined in Refs. [26,55]. We find 

i/ = 0.686(4), 7^ = 0.026(3), 7 = 1.355(8), (3.1) 

where the errors include the results of the different resummation methods. It is reassuring 
that these estimates are close to those found in Ref. [26], z/ = 0.678(10), 1] = 0.030(3), 
and 7 = 1.330(17), and also reasonably close to the Monte Carlo estimates. The small 
variation of the estimates of the critical exponents is due to the particular structure of the 
perturbative series: if they are rewritten in terms oi y = u + v and u, the resummations 
depend mostly on y* — HI, which is correctly determined by FT methods, and only shghtly 
on u* that is instead poorly known. We should also observe that the new estimate of v is 
closer to the Monte Carlo result, while the estimate of 7] is slightly worse. Therefore, the 
FT estimates do not become more accurate if more precise results for u* and v* are used. 
This is an indication that, at least for the critical exponents, the location of the fixed point 
is not the main source of error on the results. 

We also tried an alternative procedure based on an expansion of the RG functions around 
the unstable Ising fixed point u = 0, f = g*j, where [2] = 23.56(2). The analysis of the 
Ising-to-RIM RG flow reported in Ref. [28] and the discussion reported above show that it 
is convenient to introduce new variables y = u + v and z = —u. In terms of y and z, the 
RIM fixed point is located in y* ^ HI ^ 24.7(2) and z* = -3G';2 = 18.6(3), while the Ising 
fixed point is at yi — gj, z — 0. Then, we write y — Qi + Sy, obtaining for any RG function 

f{y,z) = Y.cM5y'z^, (3.2) 
Cijigi) ^Y.fijk9i- (3.3) 

k 

The value of f{y, z) at the RIM fixed point is then obtained as follows. First, we compute 
the coefficients Cij{gf) at the Ising fixed point, by using the conformal-mapping method and 
exploiting the known large-order behavior of the expansion of Cij{gi) that is determined by 
the Ising fixed point. Then, we evaluate the double series at Sy* — H^ — g} — 1.14(20) 
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TABLE X. Estimates of the coefficients Cij = (167r/3)*(67r)^Cij(c/|), cf. Eq. (3.3), for the 
expansions of rj, v, l/v, 7, and I/7 around the Ising fixed point. 







I' 


l/i/ 

/ 


7 


1/7 


Cio 


0.050(6) 


0.105(3) 


-0.273(5) 


0.181(3) 

V / 


-0.1193(8) 


C20 


0.032(8) 


0.016(4) 


0.00(1) 


0.014(5) 


0.008(6) 


C30 


0.011(6) 


0.004(5) 


0.00(3) 


0.00(1) 


0.00(1) 


Coi 





0.0500(4) 


-0.127(1) 

V / 


0.0987(6) 


-0.0646(1) 


C02 


-0.0062(2) 


0.0056(9) 


-0.003(3) 


0.015(2) 


-0.005(2) 


Cos 


0.0010(2) 


-0.003(1) 


0.017(3) 


-0.007(2) 


0.06(1) 


C04 


0.0001(5) 


0.00(1) 


0.00(1) 


0.00(2) 


0.001(1) 


Cll 





0.017(1) 


-0.003(1) 


0.032(1) 


-0.0024(2) 


C12 


-0.0018(4) 


-0.003(2) 


0.018(6) 


-0.006(3) 


0.009(3) 


Cl3 


0.0007(4) 


-0.004(3) 


0.00(1) 


-0.009(7) 


0.005(5) 


C2I 


u 


0.0033(6) 


0.0006(6) 


0.005(1) 


-0.0004(2) 


C22 


— u.uuuy ^^41 J 


-0.003(2) 


0.007(1) 


-0.005(4) 


0.003(1) 




TABLE XL Results 


obtained by 


using the expansion (3.2) for various 


truncations 





: iVlax[? + j\. ihe first 


error is due to the uncertainty on the values of the coefficients Cij{gJ), 


■ 


second one is due to (" 


ho uuc(Mlaiii(A- c 


)u tJic local ion of (lie 


Rni fixed |)oiui. 







rj-rji 


U — Uj 


1/u - l/iyj 


7-7/ 


1/7-1/7/ 


1 


0.0033(4+6) 


0.056(0+2) 


-0.145(0+5) 


0.110(1+4) 


-0.072(1+2) 


2 


-0.0025(4+8) 


0.063(1+2) 


-0.147(3+5) 


0.126(2+5) 


-0.076(1+3) 


3 


-0.0017(5+8) 


0.060(1+2) 


-0.135(3+5) 


0.120(2+4) 


-0.070(2+2) 


4 


-0.0016(6+8) 


0.062(10+3) 


-0.137(20+5) 


0.123(20+4) 


-0.070(10+2) 



and z* = — 3G22 = 18.6(3). Here, we are neglecting the fact that the RG functions are 
nonanalytic at the Ising fixed point [38-41]. Note that f{gj,0) is the value of the same 
quantity for the Ising universality class, so that the expansion (3.2) provides the differences 
between RIM and Ising critical exponents, i.e. A/ = /rim — /ising, which are expected to be 
rather small. Of course, this expansion is, at most, asymptotic. But one may hope that the 
RIM and Ising fixed points are sufficiently close, so that the first few terms of the expansion 
around the Ising fixed point allow us to obtain more accurate estimates of the exponents of 
the RIM. In Table X we report the estimates of the first coefficients Cij{g}) ior r] — r]i, v — vi, 
7 — 7/, 1/u — l/uj, and I/7 — 1/7/- The results for the critical exponents are reported in 
Table XI as a function of the order o = Max[i + j] of the expansion. We thus obtain the 
following estimates: r] - r]j = -0.0017(13), u - = 0.060(5), 1/u - l/ui = -0.135(15), 
7 — 7/ = 0.12(1), 1/7 — 1/77 = —0.07(1). The estimate we quote corresponds to o = 3, while 
the error is such to include the results with o — 2 and o — A. Then, using the estimates [9] 
777 = 0.0364(5), ui = 0.6301(4), 77 = 1.2372(5), we finally obtain 

i/ = 0.690(5), 77 = 0.035(2), 7 = 1.357(10), (3.4) 

which are again in substantial agreement with the MC results. 
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The expansion around the Ising fixed point can also be performed along the Ising-to-RIM 
RG trajectory [28], which is obtained as the limit Uo 0^ of the RG trajectories in the 
u, V plane. At least in principle, this expansion is expected to be better behaved than the 
previous one, since RG functions should be analytic near the Ising fixed point only along 
this trajectory, up to the random fixed point where nonanalyticites are again present [28] . 

An effective parametrization of the curve is given by the first few terms of its expansion 
around z — 0, which is given by 

y-yi^T(z)^C2Z^ + CsZ^ + ... (3.5) 

where [28] C2 = 0.0033(1) and Cg = 1(2) x 10"^ The fact that y - yi is of order z'^ is the 
main reason why we introduced the variable y and is due to the identity [28] 



dp. 



du 



n du 

u=0 



d/3. 



n dv 



0. (3.6) 



u=0 



Substituting the expansion (3.5) into the double expansion (3.2), we obtain an expansion in 
powers of z 

f{y,z) = f{gi + T{z),z) = Y.^i{gi)z\ (3.7) 

i 

which should be then evaluated at gi = g} and z = z* . The values of the coefficients 
ei{g}) at the Ising fixed point have been computed by using a conformal mapping and a 
Borel transform. The results for i — 1,2,3 are reported in Table XII. The estimates of 
the difference between the critical exponents of the Ising and RIM universality classes are 
r]-r]j = -0.0020(18), u - uj = 0.060(5), l/u - l/uj = -0.136(20), 7 - 7/ = 0.119(10), 
and 1/7 — 1/77 = —0.070(10). These results are obtained by truncating the expansion (3.7) 
to third order, while the error is the sum of the uncertainty due to the resummation, due 
to the truncation of the series (the difference between the second-order and the third-order 
result), and due to the uncertainty on z* (we used the Monte Carlo result). Then, by using 
the estimates [9] tji = 0.0364(5), 1/7 = 0.6301(4), 7/ = 1.2372(5), we finally obtain 

i/ = 0.690(8), 77 = 0.0345(20), 7 = 1.355(10), (3.8) 

that do not differ significantly from the estimates (3.4). 

Note that the estimate of rj obtained by using the expansion around the Ising fixed 
point is now in perfect agreement with the Monte Carlo result, at variance with the direct 
estimate (3.1). The estimate of u is also in substantial agreement with the numerical estimate 
u = 0.683(3). Therefore, the expansion around the Ising fixed point appears to be a useful 
alternative method to compute the critical properties of the RIM. 
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TABLE XII. Estimates of the coefficients Cj = (67r)*ei(5|), of. Eq. (3.7), for the expansions of 
■q, v, l/v, 7, and I/7 around the Ising fixed point. 



i 


V 


V 


\lv 


7 


1/7 


1 





0.0500(6) 


-0.1278(4) 


0.0987(6) 


-0.06462(7) 


2 


-0.0028(4) 


0.013(1) 


-0.022(3) 


0.027(2) 


-0.013(2) 


3 


0.0008(4) 


-0.002(1) 


0.012(2) 


-0.005(2) 


0.0065(8) 



APPENDIX A: NOTATIONS 



We consider the Hamiltonian (1.1) with J = 1 on a finite lattice I? with periodic 
boundary conditions. Given a quantity O depending on the spins {s} and on the random 
variable {p} we define the sample average at fixed distribution {p} 

((^>(/3, {P}) ^ E Oe-^''^^^^\ (Al) 

where Z({p}) is the sample partition function. Of course, we are interested in averaging 
over the random dilution and thus we consider disorder-averaged quantities 

(A2) 

where 

{dp\ = nM(A - 1) + (1 - x)b{p,)\. (A3) 

i 

We define the two-point correlation function and the susceptibility x(/3, V) 



G{x; (3, L) = (poso PxSx) (A4) 
x{P,L)^J2G{x;/3,L) (A5) 

X 

We also consider the second-moment correlation length ^. In infinite volume it is defined as 

eM^-^^^i:\x\'Goo{x;P). (A6) 
The finite- volume generalization is by no means unique. We use 

eiP, L) ^ GiO;P,L) -Giq^,,,,(3,L) ^ ^^^^ 

where gmin = (27r/L, 0, 0), g = 2 sing/2, and G{q] (3, L) is the Fourier transform of G{x] (3, L). 
This finite-volume definition has the correct infinite-volume limit and shows a fast conver- 
gence as L — >• cxD [49,50]. 

We also define the energy £'(/3, L) and the specific heat C(/3, L): 
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E{P,L) = 3G{e;P,L), 

(A8, 

where e = (1,0,0). We also compute higher-order couphngs. Setting 

l^k = { {^PiSi mk,k2...kn = I^kll^k2-I^K, (A9) 

i 

we define the connected n-point susceptibihties Xn averaged over random dilution by 

Vx2{P,L) = Vx{P,L)=m2, (AlO) 

VXAiP,L) = m4 - 377122, 

VxeiP, L) =mQ- 15m^2 + 30m222, 
where ^ = is the volume. Moreover, we also define 

Vx22i(3,L) = m22 - ml, 

^X42(/5, L) = 77142 - 77747772 - 3777222 + 377722^772, 

VX222{l3, L) = 777222 " 377722^772 + 27772. (AH) 

Then, we define the four-point couplings 



G22W,L) = 

and the six-point universal ratios 



X22 



exV 

H4P, L) = G4 + 3^22, (A12) 



r6(AL) = 10-^, 

xi 

X42X2 



G^2{f3,L) = A- 



X4X22 ' 



C222W,L)^6-^^. (A13) 
X22 

We will be interested in the critical value of these quantities. If S{(3, L) is any of them, we 
compute its fixed-point value (note that the order of the limits cannot be interchanged) 

S* = lim lim S{p,L). (A14) 
Finally, we define the Binder parameters 

77l„ 



n/2' 
7772' 



U22{P,L)^^, (A15) 

7772 



and the corresponding critical-point values 

U* = lim 1 

Note that the order of limits is reversed with respect to Eq. (A14). 



U* = lim lim U{/3,L). (A16) 



29 



REFERENCES 



[1] D. P. Belanger, Brazilian J. Phys. 30, 682 (2000) [cond-mat/0009029]. 

[2] M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, Phys. Rev. E 65, 066127 (2002); 

Phys. Rev. E 60, 3526 (1999). 
[3] S. Fishman and A. Aharony, J. Phys. C 12, L729 (1979). 

[4] J. L. Cardy, Phys. Rev. B 29, 505 (1984). 

[5] D. P. Belanger, in Spin Glasses and Random Fields, edited by A. P. Young (World 

Scientific, Singapore, 1998), p. 251 [cond-mat/9706042]. 
[6] The corresponding crossover exponent has been determined experimentally, finding, e.g., 

4> = 1.42(3) in Fea;Zni_2.F2 [I. B. Ferreira, A. R. King, and V. Jaccarino, Phys. Rev. B 

43, 10797 (1991)]. The theoretical estimate [P. Calabrese, A. Pelissetto, and E. Vicari, 

cond-mat/0305041] = 1.42(2) is in good agreement. 
[7] A. Aharony, in Phase Transitions and Critical Phenomena, edited by C. Domb and 

J. Lebowitz (Academic Press, New York, 1976), Vol. 6, p. 357. 
[8] R. B. Stinchcombe, in Phase Transitions and Critical Phenomena, edited by C. Domb 

and J. Lebowitz (Academic Press, New York, 1983), Vol. 7, p. 152. 
[9] A. Pehssetto and E. Vicari, Phys. Rep. 368, 549 (2002) [cond-mat/0012164]. 
[10] R. Folk, Yu. Holovatch, and T. Yavors'kn, Uspekhi Fiz. Nauk 173, 175 (2003) [Physics 

Uspekhi 46, 175 (2003)] 
[11] W. Selkc, L. N. Shchur, and A. L. Talapov, in Annual Reviews of Computational Physics, 

Vol. 1, edited by D. Stauffer (World Scientific, Singapore, 1995). 
[12] S. Wiseman and E. Domany, Phys. Rev. Lett. 81, 22 (1998); Phys. Rev. E 58, 2938 

(1998). 

[13] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. Muiioz Sudupe, G. Parisi, and 

J. J. Ruiz-Lorenzo, Phys. Rev. B 58, 2740 (1998). 
[14] K. Hukushima, J. Phys. Soc. Jpn. 69, 631 (2000). 

[15] M. I. Marques, J. A. Gonzalo, and J. Ifiiguez, Phys. Rev. E 62, 191 (2000). 

[16] W. C. Barber and D. P. Belanger, J. Magn. Magn. Mater. 226-230, 545 (2001). 

[17] P. E. Berche, C. Chatelain, B. Berche, and W. Janke, Comp. Phys. Comm. 147, 427 

(2002). 

[18] G. Grinstcin and A. Luther, Phys. Rev. B 13, 1329 (1976). 

[19] V. S. Dotsenko, B. Harris, D. Sherrington, and R. Stinchcombe, J. Phys. A 28, 3093 

(1995); V. S. Dotsenko and D. E. Feldman, J. Phys. A 28, 5183 (1995); V. S. Dotsenko, 

J. Phys. A 32, 2949 (1999). 
[20] V. V. Prudnikov, P. V. Prudnikov, and A. A. Fcdorenko, Phys. Rev. B 63, 184201 

(2001); J. Phys. A 34, L145 (2001); A. A. Fedorenko, J. Phys. A 36, 1239 (2003). 
[21] A. J. Bray, T. McCarthy, M. A. Moore, J. D. Reger, and A. P. Young, Phys. Rev. B 

36, 2212 (1987); A. J. McKane, Phys. Rev. B 49, 12003 (1994). 
[22] G. Alvarez, V. Martin-Mayor, and J. J. Ruiz- Lorenzo, J. Phys. A 33, 841 (2000). 
[23] B. N. Shalacv, S. A. Antonenko, and A. I. Sokolov, Phys. Lett. A 230, 105 (1997). 
[24] R. Folk, Yu. Holovatch, and T. Yavors'kn, Phys. Rev. B 61, 15114 (2000). 
[25] V. Blavats'ka and Yu. Holovatch, J. Mol. Liquids 105, 221 (2003). 
[26] A. Pehssetto and E. Vicari, Phys. Rev. B 62, 6393 (2000). 

[27] Perturbative field theory predicts u = 0.25(10) (Ref. [26]), u = 0.372(5) (Ref. [53]). 



30 



Moreover, it is possible to compute the next-to-leading correction exponent 002- The 

analysis of the six-loop perturbative series of Ref. [26] gives UI2 = 0.8(2) (Ref. [28]). 
[28] P. Calabrese, P. Parruccini, A. Pelissetto, and E. Vicari, cond-mat/0307699. 
[29] J.-H. Chen, M. E. Fisher, and B. G. Nickel, Phys. Rev. Lett. 48, 630 (1982). 
[30] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, and A. Munoz Sudupe, Phys. 

Lett. B 441, 330 (1998). 
[31] M. Hasenbusch, K. Pinn, and S. Vinti, Phys. Rev. B 59, 11471 (1999); M. Hasenbusch, 

J. Phys. A 32, 4851 (1999); J. Phys. A 34, 1 (2001). 
[32] M. Campostrini, M. Hasenbusch, A. Pehssetto, P. Rossi, and E. Vicari, Phys. Rev. B 

63, 214503 (2001); B 65 (2002) 144520. 
[33] S. Caracciolo, R. G. Edwards, S. J. Ferreira, A. Pelissetto, and A. D. Sokal, Phys. 

Rev. Lett. 74, 2969 (1995). 
[34] S. Caracciolo, R. G. Edwards, A. Pehssetto, and A. D. Sokal, Phys. Rev. Lett. 75, 1891 

(1995). 

[35] G. Mana, A. Pelissetto, and A. D. Sokal, Phys. Rev. D 55, 3674 (1997). 

[36] M. Liischer, P. Weisz^ and U. Wolff, Nucl. Phys. B 359, 221 (1991). 

[37] J.-K. Kim, Europhys. Lett. 28, 211 (1994); Phys. Rev. D 50, 4663 (1994); Nucl. Phys. 

B (Proc. Suppl.) 34, 702 (1994). 
[38] B. G. Nickel, in Phase Transitions, edited by M. Levy, J. C. Le Guillou, and J. Zinn- 

Justin (Plenum, New York-London, 1982). 
[39] A. Pehssetto and E. Vicari, Nucl. Phys. B 519, 626 (1998). 

[40] P. Calabrese, M. Caselle, A. Ceh, A. Pelissetto, and E. Vicari, J. Phys. A 33, 8155 
(2000). 

[41] M. Caselle, A. Pelissetto, and E. Vicari, in Fluctuating Paths and Fields, edited by W. 
Janke, A. Pelster, H.-J. Schmidt, and M. Bachmann (World Scientific, Singapore, 2001) 

[hep-th/0010228]. 

[42] P. Calabrese, M. De Prato, A. Pelissetto, and E. Vicari, cond-mat/0305434. 
[43] C. Parisi and F. Rapuano, Phys. Lett. B 157, 301 (1985). 

[44] The Parisi-Rapuano generator [43] is defined in terms of the sequence apR^„ defined 
by (ipR,n — mod (apR^„_24 + apR,n-55) 2^^). A random integer between and 2^^ — 1 
is obtained by taking rpR „ = apR „ X0RapR „_6i- We also consider the congruential 
generator r^n = mod (16807rc,„_i, 2^^ — 1). The random numbers we use are obtained 
from Tn = mod(2rc,„ + rpR,„, 2^^). 

[45] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. Muiioz Sudupe, G. Parisi, and 
J. J. Ruiz-Lorenzo, Nucl. Phys. B 512, 681 (1998). 

[46] A. Aharony and A. B. Harris, Phys. Rev. Lett. 77, 3700 (1996). 

[47] M. Palassini and S. Caracciolo, Phys. Rev. Lett. 82, 5128 (1999). 

[48] From a practical point of view, not any basis in the space spanned by the functions {z^} 
is suitable to perform the fit. As discussed in Ref. [35], care must be taken to choose 
functions which are effectively orthogonal with respect to the scalar product induced by 
the x^, otherwise the numerical solution of the minimization problem can be unstable. 
We used the functions fk{z) = yP^^°'f (2y/ymax - 1), where y = e"^/^ ?/max = e"^/^-"^'', 
P^"''') are Jacobi polynomials which are an orthonormal system in ^^^([0, l],a;"(l — xY), 
and Zmax is the maximum value of ^{(3, L)/L for the data (in our case z^^^ ~ 0.542719). 
The choice of the parameters a, 6 is not crucial and different choices would perform 



31 



analogously. 

[49] S. Caracciolo and A. Pelissetto, Phys. Rev. D 58, 105007 (1998). 

[50] With our definition of ^{(3,L) the small-2; behavior of F^{z) is determined [49] by the 
small-momentum behavior of the inverse of the two-point function of the fundamental 
fields. As in the Ising model without disorder (see V. Martin-Mayor, A. Pelissetto, and 
E. Vicari, Phys. Rev. E 66, 026112 (2002) and references therein), also in the RIM such 
corrections are extremely tiny [P. Calabrese, V. Martin-Mayor, A. Pelissetto, and E. 
Vicari, in preparation]. 

[51] D. V. Pakhnin and A. I. Sokolov, Phys. Rev. B 64, 094407 (2001); D. V. Pakhnin, 
A. I. Sokolov, and B. N. Shalaev, Pis'ma Zh. Eksp. Teor. Fiz. 75, 459 (2002) [JETP 
Lett. 75, 387 (2002)]. 

[52] Yu. Holovatch and T. Yavors'kii, J. Stat. Phys. 92, 785 (1998). 

[53] R. Folk, Yu. Holovatch, and T. Yavors'kii, Pis'ma v Zh. Eksp. Teor. Fiz. 69, 698 (1999) 

[JETP Lett. 69, 747 (1999)]. 
[54] K. B. Varnashev, Phys. Rev. B 61, 14660 (2000). 
[55] D. V. Pakhnin and A. I. Sokolov, Phys. Rev. B 61, 15130 (2000). 
[56] M. Tissier, D. Mouhanna, J. Vidal, and B. Delamotte, Phys. Rev. B 65, 140402 (2002). 
[57] J. M. Carmona, A. Pelissetto, and E. Vicari, Phys. Rev. B 61, 15136 (2000). 
[58] In Ref. [26], results are given for the rescaled couplings u and v. The relation is: u = Gttu, 

V = 167r^;/3. 



32 



